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In the recent Comment [V. M. Krasnov, Phys. Rev. B 65, 096503 (2002)], the author claims to 
have "disproved" our theoretical conclusion [S. V. Kuplevakhsky, Phys. Rev. B 60, 7496 (1999); 
ibid. B 63, 054508 (2001)] that isolated Josephson vortices in layered superconductors and stacked 
junctions are absolutely unstable in the presence of an external field. We show that this claim has 
no grounds. Moreover, by solving an appropriate boundary value problem, we obtain a complete 
classification of soliton (vortex) solutions to coupled static sine-Gordon equations. We also discuss 
the problem of vortex penetration and analyze available experimental data. 

PACS numbers: 74.50. +r, 74.80.Dm 

Unfortunately, the criticism contained in the recent Comment by_Krasnov[l] is mostly a result of a misunderstanding 
of major physical ideas and mathematical methods of my papers.Bta And wJLat is more, "the three main arguments 
against fluxons" , ascribed to me in the Comment, cannot be found in Refs.Ba: (i) my proof of the .absence of single- 
vortex solutions does not appeal to the " argument of lower free energy" ; (ii) my remark in Ref.0 on the problem 
of the existence and uniqueness concerns the difference between a finite and an infinite sets of coupled static sine- 
Gordon (SG);Q (iii) according to Refs.Htl , isolated Josephson vortices do not penetrate layered superconductors in a 
static, homogeneous external field H > .because they cannot exist inside, not vice versa. In view of this obvious 
misrepresentation of the approach of Refsrm , I have to remind the key results of these papers and the methods of 
their derivation. 

The actual main arguments of my papers are exact minimization of the Gibbs free-energy functional of the infinite 
(in the layering direction) layered superconductor in the framework of both the Lawrence-Doniach (LD) model (Ref.cl) 
and the true microscopic modelfl In particular, I resolve a nontrivial mathematical problem of the minimization with 
respect to the phases of the superconducting (S) layers. The necessity of this procedure was realized earlierJjDO.ETEl 
However, the fact that the phases and the vectot .potential are subject to an intrinsic constraint relationotl (the 
current-conservation law) was not noticed in Refs.ETu , which led to a failure. The clarification of this major issue in 
my papgfS required the use of rather sophisticated and rigorous methods of mathematical physics. As is emphasized 
in Refsrltl , the general equations derived therein encompass the whole physics of layered superconductors in parallel 
magnetic fields in the range < H < Ha .including all well-known, experimentally verified limiting cases. 

In contrast to the authors of Refs.Eni3'En3 (and the author of the Comment too), I do not make any a priori 
assumptions about the vortex structure: Vortex-plane solutions along with the Meissner solution appear as a unique 
result of exact minimization of the free-energy functional, whi ch a utomatically proves their stability and the absence 
of any other vortex solutions at H > 0. Recent experimentstiHl3 on artificial stacked junctions clearly confirm the 
existence and stability of coherent Josephson- vortex configurations that can be unambiguously identified with vortex 
planes > 0. Exceptional stability of vortex planes in the dynamic regime is confirmed by numerical simulations 
in RefsJiSLa (the authors of these papers use the terms "in-phase" and "phase-locked" solutions). n |-| 

Unfortunately, the method of exact minimization of the free-energy functional, employed in Refs.au , was not 
understood by the author of the Comment. We therefore present an alternative, simpler method in Appendix A. 

The claims of the author of the Comment that he has "disproved" our "statements" have no grounds for one 
more serious reason. In contrast to Refs.otj , the author of the Comment concentrates on the discussion of finite 
Josephson-junction stacks_and commits a typical mistake: Based on his own numerical calculations of the shape 
of the fluxon for H = 0,E_3 he makes unjustified extrapolation of the results of Refill to the case H > 0. The 
existence of single- vortex configurations at H > is assumed as " obvious" , without any mathematical proof (and 
even without an exact mathematical definition). However, the problem of the classification of vortex configurations in 
finite Josephson systems belongs to the theory of coupled SGUequations and cannot be solved by means of simplistic 
"energy arguments", numerical simulations or combinatorics. El 

In section II, we show that our conclusionscm hold for finite Josephson-junction stacks as well: Single-vortex 
configurations do not exist at H > in this case, either. However, a rigorous proof of this fact is more involved than 
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in the case of infinite layered superconductorsatl and requires different mathematical techniques .llSEJ (For this reason, 
finite Josephson-junction stacks were not considered in Refs.Ha .) Given that analytical properties of coupled-SG, 
equations have not been studied in any previous publications, we provide a necessary mathematical backgroundtatZl 
in section I. In sections III, IV, we prove the inconsistency of the rest of the arguments of the Comment. In parallel, we 
clarify misleading statements of the author of the Comment on certain important physical and mathematical issues. 

For definitiveness, we will adhere to the notation of Ref.B with a minor exception: the phase difference between the 
S-layers n and n—1 will be denoted as 4> n = Lp n — ip n ~i, with ip n being the phase. The geometry of the problem is that 
of figures 1 and 2 in Refill: The layering axis is x, with p being the period and N being the number of S-layers; the y 
axis is directed along the S-layers, with — L < y < +L being the region occupied by the system [or — oo < y < +00, if 
L = 00]. The static, homogeneous external field is applied along the z axis: H = (0, 0, H > 0). The external current 
is not considered, i.e., 1 = 0. 



I. SOME BACKGROUND INFORMATION ON ANALYTICAL PROPERTIES OF COUPLED STATIC SG 

EQUATIONS 

The difference in .mathematical description of the finite and the infinite LD models, emphasized in Ref.i , is 
not a " speculation" :tH It results from translational invariance of the barrier potential and the absence of external 
boundaries in the layering direction in the latter case. The additional symmetry of the infinite model manifests itself 
in the appearance of an unphysical degree of freedom. The elimination of this unphysical degree of freedom requires 
minimization of the Gibbs free energy with respect to the phases ip n . In contrast, minimization with respect to tp n is 
prohibited for the finite model by the physical requirement that the local field be equal to H at x = 0, x = (N — l)p. 

As an illustration, we compare coupled static SG equations for <f> n in bothrthe cases. [For r(T) -C I and H -C H C 2, 
these equations constitute solvability conditions for the Maxwell equations.!™] For N — 00, they read 

^ 2 J~TT = ~2 t( 2 + e2 ) sm( ^n _ sin 4>n+i - sin^ n _i] , (1.1) 



n = 0,±1,±2,. 



(1.2) 



Physical solutions to (O), (1.2) are subject to the requirement that the local field be equal to H at y = ±L, which 
yieldsffl 



dy 



(±L) = 2epH. 



Conditions (1.3) allow only for solutions of the type 



<t>n{y) = ~<i>n{-y) + Omod (2tt) . 



(1.3) 



(1.4) 



Note that due to the symmetry relations ( |1.4| ) conditions (1.3) alone do not specify any boundary value problem for 
(1.1), (1.2): The additive constants on the right-hand side of (1.4) remain undetermined. The determination of the 
constants 0mod(27r) that define a topological type of the solution requires imposition of boundary conditions on <f> n . 
A misunderstanding of this issue misled the author of the Comment into thinking that the solution to the boundary 
valu e prob lem for static SG equations could be "nonunique". However, the actual mathematical arbitrariness of 
(1.1), ( |l.2| ), related to the unphysical degree of freedom, lies in an infinite number of these equations: Foe. an infinite 
number of ordinary differential aauations one cannot formulate any existence and uniqueness theorems.ta The exact 
minimization procedure of Refsou (or, alternatively, of Appendix A) leads to relations 

<pn(y) = <Pn+i(y) = 4>(v) (1.5) 

that reduce (|l . l|) , (1.2) to the well-defined single.static SG equation, with the single Josephson length Aj.0 Unfortu- 
nately, in spite of "extensive theoretical studies"a of this equation, the fact that it admits an exact analytical solution, 
valid for any H and L, has not been realized in any previous publications: 



<t>{y) =tt{N v -1) + 2am {JL- 



K(k 2 ),k 



(1.6) 



2 



<ln| fcl? fc 



k 



H 



N„ = 2m, 



0,1,. 



(1.7) 



<f>(y) =TrN v +2a,m( JL-,k 



(1.8) 



H 

'in 



N„ = 2m + 1, 



0,1, 



(1.9) 



where am(u) and dn(w) = 4-am(u) are the Jakobi elliptic functions, and K (k 
kind. The superheating (penetration) field H s = H soo = (epXj) 1 determj 
Meissner state in the semiinfinite (along the layers) system < y < +oo 
specifies the number of vortex planes, with N v = for the Meissner state 
solution with N v > 1 is given by 



Hi -i 



m <h < h 



where -ff/v„ is determined by the implicit equation 

= (N v -t 

A/ 



1) 



Jk 

H, 



-K 



N v 



is the elliptic integral of the first 

s the upper bound of thepexistence of the 
The topological numbeiE3 N v = 0, 1, . . . 
The range of the existence of the soliton 



(1.10) 



(1.11) 



with Hq = H s l > H s being the superheating (penetr atio n) field for L < oo. The lower bound in ( 
appearance of the soliton) is determined by applying (1.3) to the solution of the boundary value pro 
topological conditions 

<f>(-L)=0, 0(0) =irN v , 



1.10) (the 



Diem with 



(1.12) 



whic h corresponds to the requirement that the density of the Josephson energy [and, ess entially, the free energy: see 
( |A2| )] be a minimum at the boundaries y = ±L. [Owing to the symmetry relations ( |l.4| ), it is convenient to impose 
a topological condition at y = instead of +L. The role of topological boundary conditions is exhaustively 
discussed in the literature on soliton physics.E3] The topological solution becomes unstable, when the density of the 
Josephson energy is a maximum at y = ±L, i.e., 



-L) 



0(0) = ttN v 



(1.13) 



which determines, by (1.3), the upper bound in ( 1 . lCf ) . [For the Meissner solution, an analogous condition, determining 
H s l, is obtai ned from (1.13) by setting N v = 0.] In the intermediate field range, the appropriate boundary conditions 
follow from ( |l.3| ) and continuous dependence on H: 



-f (-L) = 2e P H, 
dy 



4>(0) = ttN v . 



(1.14) 



[For the Meissner solution, we should set iV„ = in (L1J).] 



Under corresponding redefinition of the parameters A,/, H s , the solution (1.6)- (1.11) describes the single junction 
(N = 2) and the double-junction stack-(iVn= 3, see section II), which shows that vortex planes are a direct general- 
ization of ordinary Josephson vort ic es. (In the single-junction case, N v should be interpreted as the number of 
ordinary vortices.) Equations ( |l.6| )-( |Ll~l| ) cover all well-known limiting case s.E3ti m 

An important feature of the exact solution is an overlap of the regions ( 1.1 0| ) with different topological numbers 
N v . For the single junction, this overlap (not the allegedly " ncpuniqueness of the solution") was first established by 
a numerical analysisEd It is also qualitatively discussed in Ref.E3 . Mathematically, the overlap is related to the fact 
that the solution with N v cannot be continuously transformed into the solution with N v + 1 by changing H . As shown 
in section II, the overlap i s ty pical of vortex-plane solutions. 

For N < oo, equations (1.1) still hold, however, condition (fL2|) should be replaced by the conditions 



n = 1, . 



,JV-1, 



>N 



(2<AT<oo). 



(1.15) 
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Relations ( |l.3[ ), ( |l.4| ) also hold, however, relations (|ll^) should be replace by 



4>n = 0jV- 



(1.16) 



obtained from the boundary conditions on the local field at x = 0, x = (N — l)pJl3'EZl In contrast to the case N = oo, 
we can now take advantage of powerful analytical methods of the theory of ordinary differ entia l equa tions. Using 
Picard's uniqueness and existence theorem,E£l we establish the central analytical property of (LI), ( 1.15 ): 

Lemma. The initial value problem for (1.1), (1-15) with arbitrary initial conditions (f) n {yo) — ot n , -j^-iyo) — fin 
has a unique solution in the whole interval — oo < y < +oo. This solution has continuous derivatives with respect to 
y of arbitrary order and continuously depends on the initial data. 

This fundamental e xist en ce an d uniqueness Lemma will allow us to obtain a complete classification of all soliton 
(vortex) solutions to (1.1), (1.15). 



II. THE CLASSIFICATION OF SOLITON (VORTEX) SOLUTIONS TO STATIC COUPLED SG 

EQUATIONS FOR N < oo 



By analogy with (1.6)-( 1.11 ), we should find solutions to (LI), ( 1.15|) subject to topological boundary conditions 

<t> n {-L) = 0, <f> n (0) = Omod (tt) , (2.1) 

with all possible nonzero sets of the constants 0mod(7r), compatible with the requirement ( jl.3| ) that all (-L) be 
equal to each other. This task can be reformulated in terms of the initial value problem at y — —L: 

= 0, ^ {-L) = 2epH > 0, 
ay 



(2.2) 



where H is so far unknown. According to the Lemma of section I, the sol ution to (2.2), satisfying topological conditions 
at y = 0, will represent in [— L, L] the unique type of the solution to (2.1), compatible with the boundary conditions on 



-J^-. Using continuous dependence of this solution on H, we can establish a full range of its existence. This program 
is carried out to a full extent in Ref.0 . Here we present the main results. As could be expected, the only possible 
soliton solutions for H > 0, L < oo turn out to be of the vortex-plane type. In addition to the symmetry relations 
(1.4), (1.16), they satisfy the conditions 



vr < K (-L) < 0, ^(0) = 7rJV O) -p (y) > 0, n = 1, . 

ay 



.,N-1. 



(2.3) 



The range of the existence of the solution with N v > 1 is given by ( 1.10 ), where the superheating (penetration) field 



H s is now an explicit analytical function of N [H s (N < oo) > H s (N = oo) 
by the solution to the boundary value problem 

0„(-L) = -7T, 4> n (0) = irN v , 71 = 1, 



,N-1. 



- 1 ! 16 1 



and Hm is determined 



(2.4) 



Using analogous considerations, it is straightforward to establish the absence of any soliton solutions (including-the. 
vortex planes ) for H = 0, L < oo. (The latter result is indirectly confirmed by the numerical simulations of RefsJi3'E3 
for H = 0, performed in a finite interval on the y axis: Appropriate topological boundary conditions arc not fulfilled 
in the figures presented therein.) 

As a simple illustration of the above mathematical techniques, we present an explicit proof of the absence of single- 
vortex solutions to a pair of coupled static SG equations (N = 3) for H > 0. In this case, the boundary value problem 
for a single vortex, positioned at n = 1, is formulated as follows: 



, 2 (-L) = 0, 0! (0) = TT, 2 (0) = 0. 



(2.5) 



However, this boundary value problem has no solution, meeting the requirement (— L) = =^ (— L), for any H > 



dy 



Indeed, the initial value problem (2.2) for N — 3 has an exact analytical solution on the whole axis — oo < y < +oo, 
valid for any H > 0: 



</>l (y) = h(y) = 4>(y) =-n + 2am 



L 



kX 



K (k 



(2.6) 
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\ji = 



1+lepX.nH 



where Aji is the unique Josephson length for N = 3. (Contrary to a remark in the Comment, both the Meissner 

cterized by a distributi on o f Josephson lengths, whose number 
According to the Lemma, expression (2J3) is the unique solution, compatible 



solution and the vortex planes foiJl£-.N < oo are characterized by a distributi on o f Josephson lengths, whose number 
is equal to the integer part of ■y.Ej'EZl) 



with the boundary conditions at y 



-L for any H > 0. As anticipated by (2J3), it satisfies the condition > 
on the whole axis y, whereas (1.3) and (|]^) imply (0) < 0. Note that although the equality <j>\ = <p2 could be 

= 4>2 for H > was 



expected by (1.16), these symmetry relations were not employed in our proof. (The fact that <f>i 
verified experimentally on artificial double-junction stacks.til) 

Although vortex planes are unique soliton solutions for H > 0, L < oo, the situation changes drastically at H = 0, 
L = oo, when there is no need to bother about boundary conditions on -J^-. The imposition of topological boundary 
conditions 



<t>n{y) 



0, <t> n {y) = mod(7r) 



automatically ensures the fulfillment of the asymptotic boundary conditions 



dy 



(y) 



0, 



(2.7) 



(2.8) 



by virtue of Eqs. (1.1), (1-15) themselves and some ele ment ar y the orems of mathematical analysis. As shown in 
RcfsEsO , apart from the vortex-plane solution, equations (hi), ( 1.15 ) admit for H = 0, L = oo a variety of topological 
configurations: single-vortex, vortex-vortex and vortex-antivortex solutions. However, all these solutions are subject 
to the requirement that the constants 0mod(7r) be chosen from the set 0, ±7r, i.e., each 4> n can "accommodate" for 
H = 0, L = oo no more than one vortex or antivortex. (This is a typical result of thelhcQoz. of nonlinear ordinary 
differential equations that could not be obtained by means of numerical simulationsJi3'E3li3 restricted to a finite 
interval on the y axis.) In particular, the solution representing a single vortex, positioned at n = I, is specified by the 
conditions 



i(y) 



o. 



-p- (y) 

ay y- 



0, 



K0) 



dy 



(0) > 0, 



o. 



dy v- 



0, <£ n (0) = 0, 



^ (0) < 0, for n ? I, 
dy 



(2.9) 



plus a certain existence condition on all 
of ' 



(0)00 



The main distinctive feature of this solution, the negative sign 
. In this context, the remarkEI that the numerical 



X^-LO) in (P-Sl), is clearly reproduced by figures 1, 2 in Ref.b . m tms context, tnc remark 
solutionEi'^agrees well with Bulaevskiii and Clem-CoffeyS solutions" is inappropriate: Instead of exact Eqs. (|1 . 1|) , 
referencesErli3 employ a mathematically ill-defined "continuum-limit approximation" that completely neglects this 
feature. (As shown by Farid,0 "solutions" of this type are, in fact, a mathematical fiction.) 

Summarizing, the actual domain of the existence of the single- vortex solution is restricted to H = 0, L = oo. 
Concerning figure 1 in the Comment, based on the numerical calculationsEj for H = 0, the ratio of the self-energy 
of the vortex plane to that of the single vortex at H = is unimportant, in light of the absence of the latter at 
H > 0. For the same reason, the expression!!] "_ff c i = 4eTrE S i ng /<&o" is physically senseless. The actual lower critical 
field H c i is determined by the condition that the Gibbs free energy of the state with a single vortex plane be equal 
to that of the Meissner state. In the two exactly solvable cases, N = oo and N = 3 with (epHg)^ 1 « L < 
is given by H c \ = ^H s , as for the single junction.ES (In general, \J H* L — H% < H c \ < H s l.) Calculations 
of the number of " quasiequilibrium fluxon modes" by means of combinatorics are also invalid, because they do not 
take into account the necessity to satisfy topological boundary conditions together with the boundary conditions on 
(Figure 6 in Ref.c3 does not provide any evidence that such conditions for isolated fluxons are fulfilled in_the 



numerical simulations, either.) The comparisoni of the soliton vortex-plane solutions with the laminar modeJE 
for type-II superconductors is a misunderstanding: Alternating superconducting and normal layers, envisaged by this 
model, have nothing to do with soliton physics and do not possess the property of topological stability. 

The statement!!! that "free energy of any isolated solution is twice the Josephson energy" is incorrect: See, e.g., 
([A2]). However, the self-energy of soliton solutions for H = 0, -L-p oo is indeed twice the Josephson energy, which has 
been established by the use of the first integral of (1.1) in Ref.E3 (not Ref.c3 , as is-,claimed in the Comment). (From 
a field-theoretical point of view, this fact is a manifestation of the virial theoremBS) 
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III. JOSEPHSON- VORTEX PENETRATION 



Isolated Josephson vortices cannot penetrate layered superconductors and stacked junctions at any H > 0, because 
they do not form any equilibrium cpnfigurationsB'u This problem cannot be circumvented by any consideration of 
penetration as a "d^Hamic process'u whose final stage should be an equilibrium flux configuration. The experiments 
on artificial stackaI3t3 have clearly shown that Josephson vortices penetrate all the junctionSpSimultaneously and 
coherently (i.e., in the form of vortex planes), in complete agreement with the scenario of Rcijdatl . 

A description in terms of static equations, accepted in RefsBtl and previous publicationsEIrEa on single Josephson 
junctions, implies that a flux relaxation time is much smaller than a characteristic time of the change of the external 
field, which corresponds to experimental conditions of thermodynamic equilibrium. As in the case of Abrikosov 
vortices in type-II superconductors ,EZI penetration occurs when vortex planes appear at the surface y = ±L (i.e., at 
H = Hn v ), which signifies the vanishing of the surface barrier. (In contrast, the static SG equations with N > 3 
do not admit solutions with a single Josephson vortex at the surface.) Mathematically, penetration constituJ 
change of the topological type of the solution and should be interpreted as a series of first-order phase transition.^ 
Discontinuities related to these phase transitions manifest themselves, in particular, in jumps of magnetization.E3 
Such jumps of magnetization were indeed observed experimentally in stacked junctions.E3 

What figure 2 in the Comment actually shows is a violation of the conditions of thermodynamic equilibrium in the 
numerical simulations. The hypothesisa of 'decomposition of a "breather"' is also irrelevant. Certainly, the single, 
(N = 2) time-dependent SG equation for H = 0, L = oo has an exact analytical solution, known as the "breather" £3 
However, apart from the fact that time-dependent SG equations have nothing to do with thermodynamic equilibrium, 
the author of the Comment does not demonstrate the existence of such a solution foriV>3,i/>0,-L<oo and does 
not explain why it should "decompose". 

IV. EXPERIMENTAL OBSERVATIONS OF JOSEPHSON- VORTEX CONFIGURATIONS: SOME 

CONCLUDING REMARKS 

There is no contradiction between our conclusion about the absence of single-vortex configurations at H > and 
the pbservationsa of non-equilibrium isolated vortices in layered high-T c superconductors at H = 0. In the case of 
Ref.Ea , isolated vortices can be pinned by structural defects, because their energy is lower and the c-axis extent is 
smaller than those of vortex planes at H — .El 

MeasurementsEj of c-axis transport properties do not provide any direct information on flux distribution in stacked 
junctions. Both the multivaluedness and the aperiodicity oi I C (H) dependence, reported in RefE3 , can be explained 
by the overlap of states with different numbers of vortex planes, discussed in section I. |— . 

The Josephson-vortex structure in artificial stacked junctions for H > 0, I = was directly observed in Rcf.liil (by 
low-temperature scanning electron microscopy) and in Ref.E3 (by polarized neutron reflection). As is emphasized in 
the introduction and sections 11,1-111, all our theoretical conclusions fully agree with these qbaesvations. In contrast, 
the concept of "isolated fluxons"lJ stands in direct contradiction to the experiments of Refs.tiHi3 . 

In conclusion, a revision of some old prejudices in the theory of Josephson systems should begin with the case 
of the single junction, (N = 2). In spite of "extensive theoretical studies''^ for almost four decades, there is still 
an erroneous beliefflj that Josephson vortices "do not exist" ia-single junctions with L <C Xj, which is refuted by 



the exact, closed-form solution ( 1 .6 )- ( 1 . 1 1 ) . [As shown in Refs.clu , exactly the ordinary Josephson vortices and the 
vortex planes account for the occurrence of the well-known Fraunhofer pattern of I C {H) for small L.] The lack of 
the understanding of analytical properties of the single SG equation is one of the reasons why the solution to the 
more complex coupled SG equations was not pohtaiaed in previous publications. Adequate mathematical techniques 
to tackle this problem are proposed in Rcfsfl'trEjO . The author of the Comment himself could profit from the 
application of these techniques to his own research. 

APPENDIX A: A NEW METHOD OF EXACT MINIMIZATION OF THE GIBBS FREE ENERGY OF 

THE INFINITE LD MODEL 



Equations (LI), (1.2) have the first integral 



C{H) - 2^ cos <f> n (y) = Goo (n, m) — — ^— , (Al) 



G 



C(H) = 2{e P XjH) 2 — + ^cos0„(L), 

n 



Goo (n, m) 



. |n— m| 



The term on the right-hand side of (Al) is the density of electromagnetic energy. For r(T) -C 1 and H -C 7f C 2, one can 
elimi nate the electromagnetic-energy term from the LD free-energy functional^ by the use of the Maxwell equations 
and (|Al|), obtaining 



H 2 C (T)V 

47T 



- 2 +r(T) 



1 + 4 {epXjHf + ^ Yl cos ( L ) 



J dycosMv) - Aep V ' LH Y 



(A2) 



Taking the variation of ( |A2| ) with respect to tf n (y), one immediately arrives at the relations (|l.5|) . (Mathematically, 
the elimination of the electromagnetic-energy ter m is equivalent to the elimination of the intrinsic constraint in Ref.ta.) 
Unfortunately, the existence of the first integral (Al) was not noticed in any previous publications on the infinite LD 
model,BH3'ETEl which partly explains the difficulties with the minimization of the Gibbs free energy with respect to tp n . 
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